Triplet Josephson effect with magnetic feedback 
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We study AC Josephson effect in a superconductor-ferromagnet heterostructure with a variable 
magnetic configuration. The system supports triplet proximity correlations whose dynamics is 
coupled to the magnetic dynamics. This feedback dramatically modifies the behavior of the junction. 
The current-phase relation becomes double-periodic at both very low and high Josephson frequencies 
U!j. At intermediate frequencies, the periodicity in tjjt may be lost. 
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Spin-dependent transport through hybrid structures 
combining ferromagnets (F) and normal metals has at- 
tracted a lot of interest in the recent years. This interest 
is motivated by the prospect of potential technological 
applications in the field of spintronics Particular at- 
tention is given to two related effects involving mutual 
influence between the electric current through a structure 
and its magnetic configuration. The first is giant magne- 
toresistance Q in which the conductance is much larger 
when different magnetic regions have their magnetic mo- 
ments aligned than when they are anti-aligned. The op- 
posite effect is the appearance of torques acting on mag- 
netic moments when an electric current flows through 
the system These non-equilibrium current-induced 
torques appear due to non-conservation of spin currents 
accompanying a flow of charge through ferromagnetic re- 
gions. They allow manipulation of the magnetic config- 
uration, including switching between the opposite direc- 
tions or steady-state precession, without application of 
magnetic fields [4j. The two effects combined promise 
important practical applications in nonvolatile memory, 
programmable logic, and microwave oscillators. 

When the multilayer is coupled to a superconductor 
(S), an additional constraint is added, viz. that the spin 
current through the superconducting part vanishes [5(. 
This modifies the non-equilibrium torques, opening the 
possibility of perpendicular alignment of magnetic mo- 
ments. A very different situation arises when a magnetic 
structure is contacted by two superconductors. In this 
case, the proximity effect may be present, leading to a 
finite Josephson current through the structure at equilib- 
rium. The torques generated by this current correspond 
to an equilibrium effective exchange interaction between 
the magnetic moments which can be controlled by the 
phase difference between the superconductors i. The 
same mechanism enables the reciprocal effect in which 
the supercurrent depends on the magnetic configuration. 

Naive considerations might suggest that the proximity 
effect should be suppressed at short distances in the pres- 
ence of ferromagnets. However, recently it was shown 
that a long-range effect can exist due to triplet super- 
conducting correlations Q. This triplet proximity ef- 
fect (TPE), and in particular, the associated Josephson 



current, depend essentially on the magnetic configura- 
tion of the system Q. Hence S/F multilayers exhibiting 
TPE are especially suitable for studying the Josephson- 
induced magnetic exchange interaction. 

By varying the relative magnetization directions of dif- 
ferent magnetic regions, one can control the supercurrent 
flowing through the structure. Then, if the magnetic con- 
figuration is allowed to respond to the Josephson-current 
induced torques, it creates feedback for the supercurrent 
and considerably modifies it. We show that its main 
signature is frequency doubling in the current-phase re- 
lation. 

Below, we consider this feedback for a Josephson junc- 
tion biased by a dc voltage. In the AC Josephson ef- 
fect, the time dependence is normally determined by the 
current-phase relation. TPE in diffusive systems usu- 
ally leads to the conventional J — J c sin <fi relation, ex- 
cept for a special magnetic configuration with mutually 
perpendicular directions where a transition between "0" 
and "7r" states occurs Q. Then the first harmonic van- 
ishes and the current is given by the second harmonic 
~ sin 2</>; however, its amplitude is relatively small. In 
general, Josephson junctions exhibiting double-periodic 
behavior, besides being interesting objects in proximity- 
effect studies, may be useful in flux-qubit design schemes 
[§| . Josephson frequency doubling was predicted in other 
types of junctions involvingunconventional superconduc- 
tors, such as s-p [loj ]. s-d-s Q, p-p and d-d junctions with 
specific misorientation angles of the order parameter 11 1 . 
It was observed in experiments involving d-wave grain- 
boundary junctions [12]. It should be stressed, however, 
that in all these cases the frequency doubling occurs at 
isolated points in the parameter space where the first 
harmonic vanishes. Moreover, the magnitude of the cur- 
rent is suppressed in comparison with the usual value 
~ A/eRn, where R n is the normal-state resistance. 

In this work, we consider the magnetic exchange in- 
teraction induced by Josephson currents in a dirty S/F 
heterostructure exhibiting TPE. We show that this inter- 
action may prefer non-collinear magnetic configurations 
and the preferred direction depends continuously on the 
superconducting phase difference. Thus, the static mag- 
netic configuration can be controlled by the applied phase 
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difference. We then consider the influence of feedback 
from the magnetic moments on the AC Josephson effect. 
The magnetic system exhibits a range of different be- 
haviors, from simple harmonic oscillations to fractional- 
frequency periodic behavior and chaotic motion. A fi- 
nite zero-frequency deviation from the equilibrium con- 
figuration may appear, allowing control of the direction 
of the average magnetization also by an applied volt- 
age. The magnetic feedback complicates the behavior 
of the current in the time domain, making it generally 
impossible to express it in terms of a current-phase re- 
lation. On the other hand, we find that both in the 
low- and high frequency limit such a relation becomes 
meaningful, with the current exhibiting a double-phase 
dependence, J ~ sin20(i) or J ~ cos20(i). The criti- 
cal current in the low-frequency regime is of the order of 
the value Exh/sRn, characteristic for diffusive systems. 
The unusual cosine dependence of the Josephson current 
appears when Gilbert damping is important in the mag- 
netic dynamics, breaking the time-reversal symmetry. At 
high frequencies, the magnetization cannot effectively fol- 
low the phase variation, leading to a ~ l/u> 2 suppression 
of the effective Josephson coupling. At even higher fre- 
quencies, the damping is dominant, and the frequency 
dependence becomes ~ 1/lu. The presence of damping 
is expressed in the appearance of a dc component of the 
current leading to a finite resistance. 

The system. We consider an S/F heterostructure 
described in Fig. [1] which is a minimal discrete setup ex- 
hibiting the triplet proximity effect. Two magnetic re- 
gions 1 and 3 are adjacent to the superconducting reser- 
voirs that induce proximity mini-gaps Ai ( 3 in them. Be- 
tween these regions there is an additional magnetic re- 
gion 2 whose length is much larger than fa and where 
triplet superconducting correlations are induced. This 
region is assumed to be weakly polarized (metallic), so 
that both spin directions are present at the Fermi sur- 
face. The magnetic regions are characterized by the ex- 
change energies hi, while the magnetization directions 
rii are specified by the angles Ox, 9% and \ as shown in 
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Fig. [TJ Assuming that the conductances of these regions 
are much higher than the conductances g™ 3 of the con- 
nectors between them, our system can be described by a 
circuit-theory model for the triplet proximity effect used 
in Ref. 8. Magnetization directions of regions 1 and 2 are 
assumed fixed, e.g. by pinning to an antifcrromagnetic 
substrate, or by geometrical shaping, with the angle be- 
tween them being 9±. On the other hand, magnetization 
ri3 is free to rotate, with region 3 separated by a normal 
spacer from region 2 in order to avoid exchange coupling 
between them. 

In accordance with the model assumptions, regions 1 
and 3 act as effective S-F reservoirs, hence their energies 
are independent of the magnetic configuration. On the 
other hand, triplet superconducting correlations extend- 
ing through region 2, are very sensitive to the magnetiza- 
tion directions. Hence the configuration-dependent part 
of the energy can be found by integrating over the den- 
sity of states (DOS) in region 2. The DOS for each spin 
direction is given by Q 



(e) = > 1- 
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where vq is 
g%\A k \smO k /(g? 



the 
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normal-state DOS, = 



g§)y/hi-\A k 

is the superconducting phase differ- 



ence, and Eth is the Thouless energy of the structure. 
Using this expression, one can see that the energy is 
given by a logarithmic integral and the main contribu- 
tion comes from e 3> Eth- In the leading order one 
obtains 



vqv 2 A cut . 
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2ai<2,3 cos 4> cos x) 



(2) 

where vi is the volume of the magnetic region 2 and 
A cut ~ min(Ai,/ii — A,-) is a cutoff energy. This ex- 
pression can be written in a form presenting explicitly 
the dependence on the orientation angles 63 and Xi 

E =p\ sin 2 6*3 + 2pip 3 sin #3 cos 4> cos x, (3) 

with P1.3 being effective exchange couplings for the mag- 
netic vector TI3. The stable configuration is achieved 
when all magnetization directions are in the same plane, 
denoted in the following as the x — z plane, and n-j is 
tilted with respect to rii by a finite angle satisfying 



sin 03 = — I cos ( 

P3 



(4) 



FIG. 1: The experimental setup 



This angle depends continuously on the applied super- 
conducting phase difference <f>, while the angle x assumes 
the values or it so that the product cos 4> cos x is nega- 
tive. In fact, there are two stable directions, given by the 
angles #3 and 7r — 63 . In what follows we will treat them 
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FIG. 2: Trajectories of the magnetization vector in the x — y 
plane for different Josephson frequencies (given in units of 
uj m )- Trajectory (b) is chaotic, while trajectory (c) has a 
finite zero-frequency component for n x . Here T is the period 
of the trajectory, and To = 2tt/wj. For comparison, the low- 
frequency trajectory lies entirely on the x axis. 



as equivalent, since they correspond to the same current. 
The energy of the stable configuration is given by 



E m in = ~Pl COS 2 ( 



(5) 



Hence allowing the magnetization direction 713 to orient 
itself along the stable direction leads to the current-phase 
relation J = J c sin 2<p. 

Low frequencies. When a small voltage V is applied 
to the structure, such that the corresponding frequency 
u>j = 2eV/h is much smaller than the characteristic fre- 
quency of the magnetic system ui m (see below), the vector 
77,3 follows the stable direction given by Eq. (j4]) , perform- 
ing slow oscillations in the x — z plane. The alternating 
Josephson current oscillates with the double frequency 



J = —Pi sm ■ 
a 



AeV 



-t , 



(6) 



while the critical current remains of the same order of 
magnitude as in the case with a fixed magnetic configu- 
ration. 

For higher Josephson frequencies, the variation of 77,3 
is no more limited to the x — z plane. Instead, the mag- 
netization performs a variety of non-harmonic motions 
whose frequency may be a multiple or a fraction of the 
driving frequency luj [Fig. [2] (a), (c), (d)]. For certain 
trajectories the time average of 63 is finite [Fig. [2] (c)], 
corresponding to a tilt of 77,3 away from the equilibrium 
in response to an applied voltage. Within some frequency 
intervals, the motion is chaotic, as shown in Fig.[2](b). In 
these intermediate regimes, the Josephson current shows 



a complicated time dependence which is generally not 
periodic in 2tt/ujj. Hence this dependence cannot be pa- 
rameterized in terms of the phase. Instead, one can speak 
of a Josephson current with a time-dependent coupling. 

High frequencies. When applied voltage is high, the 
Josephson frequency becomes much higher than the mag- 
netic frequencies. In this case the magnetic vector 77,3 
cannot effectively follow the fast oscillations of the po- 
tential, and the time-averaged potential seen by 77,3 has 
a minimum for 723 || z. The motion of 713 can be deter- 
mined by expanding 713 = z + 5n and using a linearized 
Landau-Lifshits-Gilbert (LLG) equation, 



8h = z(— 7 x H e ff + aSn) 



(7) 



where 7 is the gyromagnctic ratio, a is the effective 
damping coefficient, i? e // = —dE/dm^ is the effective 
field, and is the magnetization density of magnetic 
region 3. 

When Gilbert damping is negligible, the trajectory of 
77,3 has a very low aspect ratio, so that the motion is 
almost completely confined to the y axis. It is given by 



5n, 



j 2 piplh 2 
e 2 V 2 m 2 

eVrriz 



2eVt 
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(8) 



Thus at high frequencies, 77,3 precesses in phase with the 
voltage pumping. This leads to an increase in the Joseph- 
son energy, and, correspondingly, a negative Josephson 
current, 



J 



2^ ( ^Ps ), in AeVt 
e \ Vm^ I Ti 



(9) 



Hence in the high-frequency regime the system shows not 
only frequency doubling, but also an effective 7r-junction 
behavior. The magnitude of the current is suppressed as 
~ V~ 2 as shown in Fig. [3] 

The neglect of damping is justified as long as auj <C 
w m = jp 2 /m 3 . When the voltage is high enough, this 
condition is not satisfied anymore, and the dissipation 
starts to be important. As the Josephson frequency be- 
comes so large that the opposite inequality holds, the 
motion of TI3 is determined by the driving against the 
damping force, 



. IPWih „ . 2eVt 

on = — — (—ax + y) sm — - — 

eVm 3 " n 

Then the Josephson current is given by 

4eVt\ 



j = 2a 1V \ V \ i 



cos ■ 



h J 



(10) 



(11) 



Note the unusual cosine dependence on the phase. It oc- 
curs since the time-reversal symmetry is broken by the 
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FIG. 3: The absolute value of the Josephson current harmon- 
ics proportional to sin 2u> jt (asterisks) and to 1 — cos 2uijt 
(dots). Solid lines are fits ~ 1/uj and 1/lu 2 . The data 
points are obtained from numerical integration of the full 
(non-linear) LLG equation. 

dissipation in this regime. Due to the same reason, a 
zero-frequency component of the current appears, signi- 
fying the onset of a finite nonlinear resistance across the 
structure. Since this regime is governed by the damp- 
ing, the current amplitude is proportional to a, while 
the suppression ~ 1/V is weaker in this regime (Fig. [3]). 

To estimate the magnetic dynamics frequency u> m , we 
use typical values Eth ~ 1 meV, v Q ~ l/(eV/atom), 
and 7713 ~ 1/xs/atom, where \xb is the Bohr magneton. 
Then w m ~ GHz, where ^2,3 are the volumes of 

the corresponding magnetic regions. As this frequency 
is quite low, observation of the high-frequency regimes 
should present no difficulty. On the other hand, the low- 
frequency AC regime would require extremely low volt- 
ages, below 1 /iV. A reasonable alternative would be in- 
corporating the structure in a superconducting loop and 
measuring the Josephson current as a function of the ap- 
plied flux. 

Applicability of our model requires that any mag- 
netic anisotropy of part 3 should be smaller than the 
proximity-induced energy, Eq. With the above val- 
ues of the parameters it is of the order of 10 4 x u 2 J/m 3 , 
so one should choose materials with low value of the 
cristalline anisotropy, such as permalloy. Finally, we em- 
phasize that the properties discussed above are specific 
for metallic systems. In half-metals, the behavior will 
be very different. Thus, in the low- frequency regime n.3 
precesses around n 2 at a constant angle #3, while the 
Josephson current vanishes. 

Conclusions. We have considered the AC Joseph- 
son effect in a S/F/S structure with magnetic dynam- 
ics coupled to the dynamics of superconducting correla- 
tions. The magnetic configuration in the structure was 
assumed to be non-uniform so that the structure ex- 



hibits a triplet proximity effect. Variation of the mag- 
netic configuration is shown to essentially modify the 
current behavior that can be observed in the appearance 
of fractional Shapiro steps. Thus measurement of the 
Josephson current would provide information about the 
coupling and self-consistent feedback dynamics between 
the superconducting and magnetic degrees of freedom. 
The coupling also allows to control the magnetization 
direction by means of applied voltage or superconduct- 
ing phase. In the low-frequency limit, the magnetization 
follows the immediate potential minimum, leading to a 
~ sin 20 current-phase relation. The critical current has 
the same order of magnitude Exh/sRn as that due to the 
usual singlet proximity effect in dirty structures. In the 
high-frequency regime, as long as the damping is not im- 
portant, the Josephson current is negative, correspond- 
ing to a 7r-junction behavior. It is suppressed by a factor 
~ (cu m /u!j) 2 relative to the low-frequency regime. At 
even higher frequencies, Gilbert damping starts playing 
the major role in the dynamics. Then the time-reversal 
symmetry is broken and the current-phase relation takes 
an unusual cosine form. In addition, a DC component 
of the current appears, manifesting itself in a finite resis- 
tance. The current suppression becomes weaker in this 
regime. 
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